The Phylogeny of Bacillus cereus sensu lato.

The three main species of the Bacillus cereus sensu lato, B. cereus, B. thuringiensis, and B. anthracis, were recognized and established by the early 1900 s because they each exhibited distinct phenotypic traits. B. thuringiensis isolates and their parasporal crystal proteins have long been established as a natural pesticide and insect pathogen. B. anthracis, the etiological agent for anthrax, was used by Robert Koch in the 19th century as a model to develop the germ theory of disease, and B. cereus, a common soil organism, is also an occasional opportunistic pathogen of humans. In addition to these three historical species designations, are three less-recognized and -understood species: B. mycoides, B. weihenstephanensis, and B. pseudomycoides. All of these "species" combined comprise the Bacillus cereus sensu lato group. Despite these apparently clear phenotypic definitions, early molecular approaches to separate the first three by various DNA hybridization and 16S/23S ribosomal sequence analyses led to some "confusion" because there were limited differences to differentiate between these species. These and other results have led to frequent suggestions that a taxonomic change was warranted to reclassify this group to a single species. But the pathogenic properties of B. anthracis and the biopesticide applications of B. thuringiensis appear to "have outweighed pure taxonomic considerations" and the separate species categories are still being maintained. B. cereus sensu lato represents a classic example of a now common bacterial species taxonomic quandary.

Pangenome Analysis of Burkholderia pseudomallei: Genome Evolution Preserves Gene Order despite High Recombination Rates.

The pangenomic diversity in Burkholderia pseudomallei is high, with approximately 5.8% of the genome consisting of genomic islands. Genomic islands are known hotspots for recombination driven primarily by site-specific recombination associated with tRNAs. However, recombination rates in other portions of the genome are also high, a feature we expected to disrupt gene order. We analyzed the pangenome of 37 isolates of B. pseudomallei and demonstrate that the pangenome is 'open', with approximately 136 new genes identified with each new genome sequenced, and that the global core genome consists of 4568±16 homologs. Genes associated with metabolism were statistically overrepresented in the core genome, and genes associated with mobile elements, disease, and motility were primarily associated with accessory portions of the pangenome. The frequency distribution of genes present in between 1 and 37 of the genomes analyzed matches well with a model of genome evolution in which 96% of the genome has very low recombination rates but 4% of the genome recombines readily. Using homologous genes among pairs of genomes, we found that gene order was highly conserved among strains, despite the high recombination rates previously observed. High rates of gene transfer and recombination are incompatible with retaining gene order unless these processes are either highly localized to specific sites within the genome, or are characterized by symmetrical gene gain and loss. Our results demonstrate that both processes occur: localized recombination introduces many new genes at relatively few sites, and recombination throughout the genome generates the novel multi-locus sequence types previously observed while preserving gene order.

Genome Sequence of Bacillus anthracis STI, a Sterne-Like Georgian/Soviet Vaccine Strain.

The Bacillus anthracis strain STI is a Soviet vaccine strain that lacks the pXO2 plasmid. Previous data indicate that this isolate forms a new branch within the B. anthracis sub-group originally identified as A. Br.008/009.

Phylogeography of Bacillus anthracis in the country of Georgia shows evidence of population structuring and is dissimilar to other regional genotypes.

Sequence analyses and subtyping of Bacillus anthracis strains from Georgia reveal a single distinct lineage (Aust94) that is ecologically established. Phylogeographic analysis and comparisons to a global collection reveals a clade that is mostly restricted to Georgia. Within this clade, many groups are found around the country, however at least one subclade is only found in the eastern part. This pattern suggests that dispersal into and out of Georgia has been rare and despite historical dispersion within the country, for at least for one lineage, current spread is limited.

Molecular epidemiologic investigation of an anthrax outbreak among heroin users, Europe.

In December 2009, two unusual cases of anthrax were diagnosed in heroin users in Scotland. A subsequent anthrax outbreak in heroin users emerged throughout Scotland and expanded into England and Germany, sparking concern of nefarious introduction of anthrax spores into the heroin supply. To better understand the outbreak origin, we used established genetic signatures that provided insights about strain origin. Next, we sequenced the whole genome of a representative Bacillus anthracis strain from a heroin user (Ba4599), developed Ba4599-specific single-nucleotide polymorphism assays, and genotyped all available material from other heroin users with anthrax. Of 34 case-patients with B. anthracis-positive PCR results, all shared the Ba4599 single-nucleotide polymorphism genotype. Phylogeographic analysis demonstrated that Ba4599 was closely related to strains from Turkey and not to previously identified isolates from Scotland or Afghanistan, the presumed origin of the heroin. Our results suggest accidental contamination along the drug trafficking route through a cutting agent or animal hides used to smuggle heroin into Europe.

An attenuated strain of Bacillus anthracis (CDC 684) has a large chromosomal inversion and altered growth kinetics.

BACKGROUND: An isolate originally labeled Bacillus megaterium CDC 684 was found to contain both pXO1 and pXO2, was non-hemolytic, sensitive to gamma-phage, and produced both the protective antigen and the poly-D-glutamic acid capsule. These phenotypes prompted Ezzell et al., (J. Clin. Microbiol. 28:223) to reclassify this isolate to Bacillus anthracis in 1990. RESULTS: We demonstrate that despite these B. anthracis features, the isolate is severely attenuated in a guinea pig model. This prompted whole genome sequencing and closure. The comparative analysis of CDC 684 to other sequenced B. anthracis isolates and further analysis reveals: a) CDC 684 is a close relative of a virulent strain, Vollum A0488; b) CDC 684 defines a new B. anthracis lineage (at least 51 SNPs) that includes 15 other isolates; c) the genome of CDC 684 contains a large chromosomal inversion that spans 3.3 Mbp; d) this inversion has caused a displacement of the usual spatial orientation of the origin of replication (ori) to the termination of replication (ter) from 180° in wild-type B. anthracis to 120° in CDC 684 and e) this isolate also has altered growth kinetics in liquid media. CONCLUSIONS: We propose two alternative hypotheses explaining the attenuated phenotype of this isolate. Hypothesis 1 suggests that the skewed ori/ter relationship in CDC 684 has altered its DNA replication and/or transcriptome processes resulting in altered growth kinetics and virulence capacity. Hypothesis 2 suggests that one or more of the single nucleotide polymorphisms in CDC 684 has altered the expression of a regulatory element or other genes necessary for virulence.

Phylogeographic reconstruction of a bacterial species with high levels of lateral gene transfer.

BACKGROUND: Phylogeographic reconstruction of some bacterial populations is hindered by low diversity coupled with high levels of lateral gene transfer. A comparison of recombination levels and diversity at seven housekeeping genes for eleven bacterial species, most of which are commonly cited as having high levels of lateral gene transfer shows that the relative contributions of homologous recombination versus mutation for Burkholderia pseudomallei is over two times higher than for Streptococcus pneumoniae and is thus the highest value yet reported in bacteria. Despite the potential for homologous recombination to increase diversity, B. pseudomallei exhibits a relative lack of diversity at these loci. In these situations, whole genome genotyping of orthologous shared single nucleotide polymorphism loci, discovered using next generation sequencing technologies, can provide very large data sets capable of estimating core phylogenetic relationships. We compared and searched 43 whole genome sequences of B. pseudomallei and its closest relatives for single nucleotide polymorphisms in orthologous shared regions to use in phylogenetic reconstruction. RESULTS: Bayesian phylogenetic analyses of >14,000 single nucleotide polymorphisms yielded completely resolved trees for these 43 strains with high levels of statistical support. These results enable a better understanding of a separate analysis of population differentiation among >1,700 B. pseudomallei isolates as defined by sequence data from seven housekeeping genes. We analyzed this larger data set for population structure and allele sharing that can be attributed to lateral gene transfer. Our results suggest that despite an almost panmictic population, we can detect two distinct populations of B. pseudomallei that conform to biogeographic patterns found in many plant and animal species. That is, separation along Wallace's Line, a biogeographic boundary between Southeast Asia and Australia. CONCLUSION: We describe an Australian origin for B. pseudomallei, characterized by a single introduction event into Southeast Asia during a recent glacial period, and variable levels of lateral gene transfer within populations. These patterns provide insights into mechanisms of genetic diversification in B. pseudomallei and its closest relatives, and provide a framework for integrating the traditionally separate fields of population genetics and phylogenetics for other bacterial species with high levels of lateral gene transfer.

Phylogenetic understanding of clonal populations in an era of whole genome sequencing.

Phylogenetic hypotheses using whole genome sequences have the potential for unprecedented accuracy, yet a failure to understand issues associated with discovery bias, character sampling, and strain sampling can lead to highly erroneous conclusions. For microbial pathogens, phylogenies derived from whole genome sequences are becoming more common, as large numbers of characters distributed across entire genomes can yield extremely accurate phylogenies, particularly for strictly clonal populations. The availability of whole genomes is increasing as new sequencing technologies reduce the cost and time required for genome sequencing. Until entire sample collections can be fully sequenced, harnessing the phylogenetic power from whole genome sequences in more than a small subset of fully sequenced strains requires the integration of whole genome and partial genome genotyping data. Such integration involves discovering evolutionarily stable polymorphic characters by whole genome comparisons, then determining allelic states across a wide panel of isolates using high-throughput genotyping technologies. Here, we demonstrate how such an approach using single nucleotide polymorphisms (SNPs) yields highly accurate, but biased phylogenetic reconstructions and how the accuracy of the resulting tree is compromised by incomplete taxon and character sampling. Despite recent phylogenetic work detailing the strengths and biases of integrating whole genome and partial genome genotype data, these issues are relatively new and remain poorly understood by many researchers. Here, we revisit these biases and provide strategies for maximizing phylogenetic accuracy. Although we write this review with bacterial pathogens in mind, these concepts apply to any clonally reproducing population or indeed to any evolutionarily stable marker that is inherited in a strictly clonal manner. Understanding the ways in which current and emerging technologies can be used to maximize phylogenetic knowledge is advantageous only with a complete understanding of the strengths and weaknesses of these methods.

Bacillus anthracis in China and its relationship to worldwide lineages.

BACKGROUND: The global pattern of distribution of 1033 B. anthracis isolates has previously been defined by a set of 12 conserved canonical single nucleotide polymorphisms (canSNP). These studies reinforced the presence of three major lineages and 12 sub-lineages and sub-groups of this anthrax-causing pathogen. Isolates that form the A lineage (unlike the B and C lineages) have become widely dispersed throughout the world and form the basis for the geographical disposition of "modern" anthrax. An archival collection of 191 different B. anthracis isolates from China provides a glimpse into the possible role of Chinese trade and commerce in the spread of certain sub-lineages of this pathogen. Canonical single nucleotide polymorphism (canSNP) and multiple locus VNTR analysis (MLVA) typing has been used to examine this archival collection of isolates. RESULTS: The canSNP study indicates that there are 5 different sub-lineages/sub-groups in China out of 12 previously described world-wide canSNP genotypes. Three of these canSNP genotypes were only found in the western-most province of China, Xinjiang. These genotypes were A.Br.008/009, a sub-group that is spread across most of Europe and Asia; A.Br.Aust 94, a sub-lineage that is present in Europe and India, and A.Br.Vollum, a lineage that is also present in Europe. The remaining two canSNP genotypes are spread across the whole of China and belong to sub-group A.Br.001/002 and the A.Br.Ames sub-lineage, two closely related genotypes. MLVA typing adds resolution to the isolates in each canSNP genotype and diversity indices for the A.Br.008/009 and A.Br.001/002 sub-groups suggest that these represent older and established clades in China. CONCLUSION: B. anthracis isolates were recovered from three canSNP sub-groups (A.Br.008/009, A.Br.Aust94, and A.Br.Vollum) in the western most portion of the large Chinese province of Xinjiang. The city of Kashi in this province appears to have served as a crossroads for not only trade but the movement of diseases such as anthrax along the ancient "silk road". Phylogenetic inference also suggests that the A.Br.Ames sub-lineage, first identified in the original Ames strain isolated from Jim Hogg County, TX, is descended from the A.Br.001/002 sub-group that has a major presence in most of China. These results suggest a genetic discontinuity between the younger Ames sub-lineage in Texas and the large Western North American sub-lineage spread across central Canada and the Dakotas.

Pre-Columbian origins for North American anthrax.

Disease introduction into the New World during colonial expansion is well documented and had a major impact on indigenous populations; however, few diseases have been associated with early human migrations into North America. During the late Pleistocene epoch, Asia and North America were joined by the Beringian Steppe ecosystem which allowed animals and humans to freely cross what would become a water barrier in the Holocene. Anthrax has clearly been shown to be dispersed by human commerce and trade in animal products contaminated with Bacillus anthracis spores. Humans appear to have brought B. anthracis to this area from Asia and then moved it further south as an ice-free corridor opened in central Canada approximately 13,000 ybp. In this study, we have defined the evolutionary history of Western North American (WNA) anthrax using 2,850 single nucleotide polymorphisms (SNPs) and 285 geographically diverse B. anthracis isolates. Phylogeography of the major WNA B. anthracis clone reveals ancestral populations in northern Canada with progressively derived populations to the south; the most recent ancestor of this clonal lineage is in Eurasia. Our phylogeographic patterns are consistent with B. anthracis arriving with humans via the Bering Land Bridge. This northern-origin hypothesis is highly consistent with our phylogeographic patterns and rates of SNP accumulation observed in current day B. anthracis isolates. Continent-wide dispersal of WNA B. anthracis likely required movement by later European colonizers, but the continent's first inhabitants may have seeded the initial North American populations.

Single nucleotide polymorphism typing of Bacillus anthracis from Sverdlovsk tissue.

A small number of conserved canonical single nucleotide polymorphisms (canSNP) that define major phylogenetic branches for Bacillus anthracis were used to place a Sverdlovsk patient's B. anthracis genotype into 1 of 12 subgroups. Reconstruction of the pagA gene also showed a unique SNP that defines a new lineage for B. anthracis.

Differentiation of Clostridium botulinum serotype A strains by multiple-locus variable-number tandem-repeat analysis.

Ten variable-number tandem-repeat (VNTR) regions identified within the complete genomic sequence of Clostridium botulinum strain ATCC 3502 were used to characterize 59 C. botulinum strains of the botulism neurotoxin A1 (BoNT/A1) to BoNT/A4 (BoNT/A1-A4) subtypes to determine their ability to discriminate among the serotype A strains. Two strains representing each of the C. botulinum serotypes B to G, including five bivalent strains, and two strains of the closely related species Clostridium sporogenes were also tested. Amplified fragment length polymorphism analyses revealed the genetic diversity among the serotypes and the high degree of similarity among many of the BoNT/A1 strains. The 10 VNTR markers amplified fragments within all of the serotype A strains but were less successful with strains of other serotypes. The composite multiple-locus VNTR analysis of the 59 BoNT/A1-A4 strains and 3 bivalent B strains identified 38 different genotypes. Thirty genotypes were identified among the 53 BoNT/A1 and BoNT/A1(B) strains, demonstrating discrimination below the subtype level. Contaminating DNA within crude toxin preparations of three BoNT/A subtypes (BoNT/A1 to BoNT/A3) also supported amplification of all of the VNTR regions. These markers provide clinical and forensics laboratories with a rapid, highly discriminatory tool to distinguish among C. botulinum BoNT/A1 strains for investigations of botulism outbreaks.

Real-time PCR assays of single-nucleotide polymorphisms defining the major Brucella clades.

Members of the genus Brucella are known worldwide as pathogens of wildlife and livestock and are the most common organisms of zoonotic infection in humans. In general, brucellae exhibit a range of host specificity in animals that has led to the identification of at least seven Brucella species. The genomes of the various Brucella species are highly conserved, which makes the differentiation of species highly challenging. However, we found single-nucleotide polymorphisms (SNPs) in housekeeping and other genes that differentiated the seven main Brucella species or clades and thus enabled us to develop real-time PCR assays based around these SNPs. Screening of a diverse panel of 338 diverse isolates with these assays correctly identified each isolate with its previously determined Brucella clade. Six of the seven clade-specific assays detected DNA concentrations of less than 10 fg, indicating a high level of sensitivity. This SNP-based approach places samples into a phylogenetic framework, allowing reliable comparisons to be made among the lineages of clonal bacteria and providing a solid basis for genotyping. These PCR assays provide a rapid and highly sensitive method of differentiating the major Brucella groups that will be valuable for clinical and forensic applications.

Evaluating Burkholderia pseudomallei Bip proteins as vaccines and Bip antibodies as detection agents.

Burkholderia pseudomallei is a biothreat agent and an important natural pathogen, causing melioidosis in humans and animals. A type III secretion system (TTSS-3) has been shown to be critical for virulence. Because TTSS components from other pathogens have been used successfully as diagnostic agents and as experimental vaccines, it was investigated whether this was the case for BipB, BipC and BipD, components of B. pseudomallei's TTSS-3. The sequences of BipB, BipC and BipD were found to be highly conserved among B. pseudomallei and B. mallei isolates. A collection of monoclonal antibodies (mAbs) specific for each Bip protein was obtained. Most recognized both native and denatured Bip protein. Burkholderia pseudomallei or B. mallei did not express detectable BipB or BipD under the growth conditions used. However, anti-BipD mAbs did recognize the TTSS needle structures of a Shigella strain engineered to express BipD. The authors did not find that BipB, BipC or BipD are protective antigens because vaccination of mice with any single protein did not result in protection against experimental melioidosis. Enzyme-linked immunosorbent assay (ELISA) studies showed that human melioidosis patients had antibodies to BipB and BipD. However, these ELISAs had low diagnostic accuracy in endemic regions, possibly due to previous patient exposure to B. pseudomallei.

Global genetic population structure of Bacillus anthracis.

Anthrax, caused by the bacterium Bacillus anthracis, is a disease of historical and current importance that is found throughout the world. The basis of its historical transmission is anecdotal and its true global population structure has remained largely cryptic. Seven diverse B. anthracis strains were whole-genome sequenced to identify rare single nucleotide polymorphisms (SNPs), followed by phylogenetic reconstruction of these characters onto an evolutionary model. This analysis identified SNPs that define the major clonal lineages within the species. These SNPs, in concert with 15 variable number tandem repeat (VNTR) markers, were used to subtype a collection of 1,033 B. anthracis isolates from 42 countries to create an extensive genotype data set. These analyses subdivided the isolates into three previously recognized major lineages (A, B, and C), with further subdivision into 12 clonal sub-lineages or sub-groups and, finally, 221 unique MLVA15 genotypes. This rare genomic variation was used to document the evolutionary progression of B. anthracis and to establish global patterns of diversity. Isolates in the A lineage are widely dispersed globally, whereas the B and C lineages occur on more restricted spatial scales. Molecular clock models based upon genome-wide synonymous substitutions indicate there was a massive radiation of the A lineage that occurred in the mid-Holocene (3,064-6,127 ybp). On more recent temporal scales, the global population structure of B. anthracis reflects colonial-era importation of specific genotypes from the Old World into the New World, as well as the repeated industrial importation of diverse genotypes into developed countries via spore-contaminated animal products. These findings indicate humans have played an important role in the evolution of anthrax by increasing the proliferation and dispersal of this now global disease. Finally, the value of global genotypic analysis for investigating bioterrorist-mediated outbreaks of anthrax is demonstrated.

The complete genome sequence of Bacillus thuringiensis Al Hakam.

Bacillus thuringiensis is an insect pathogen that is widely used as a biopesticide (E. Schnepf, N. Crickmore, J. Van Rie, D. Lereclus, J. Baum, J. Feitelson, D. R. Zeigler, and D. H. Dean, Microbiol. Mol. Biol. Rev. 62:775-806, 1998). Here we report the finished, annotated genome sequence of B. thuringiensis Al Hakam, which was collected in Iraq by the United Nations Special Commission (L. Radnedge, P. Agron, K. Hill, P. Jackson, L. Ticknor, P. Keim, and G. Andersen, Appl. Environ. Microbiol. 69:2755-2764, 2003).

Characterization of Bacillus cereus isolates associated with fatal pneumonias: strains are closely related to Bacillus anthracis and harbor B. anthracis virulence genes.

Bacillus cereus is ubiquitous in nature, and while most isolates appear to be harmless, some are associated with food-borne illnesses, periodontal diseases, and other more serious infections. In one such infection, B. cereus G9241 was identified as the causative agent of a severe pneumonia in a Louisiana welder in 1994. This isolate was found to harbor most of the B. anthracis virulence plasmid pXO1 (13). Here we report the characterization of two clinical and one environmental B. cereus isolate collected during an investigation of two fatal pneumonia cases in Texas metal workers. Molecular subtyping revealed that the two cases were not caused by the same strain. However, one of the three isolates was indistinguishable from B. cereus G9241. PCR analysis demonstrated that both clinical isolates contained B. anthracis pXO1 toxin genes. One clinical isolate and the environmental isolate collected from that victim's worksite contained the cap A, B, and C genes required for capsule biosynthesis in B. anthracis. Both clinical isolates expressed a capsule; however, neither was composed of poly-D-glutamic acid. Although most B. cereus isolates are not opportunistic pathogens and only a limited number cause food-borne illnesses, these results demonstrate that some B. cereus strains can cause severe and even fatal infections in patients who appear to be otherwise healthy.

Pathogenomic sequence analysis of Bacillus cereus and Bacillus thuringiensis isolates closely related to Bacillus anthracis.

Bacillus anthracis, Bacillus cereus, and Bacillus thuringiensis are closely related gram-positive, spore-forming bacteria of the B. cereus sensu lato group. While independently derived strains of B. anthracis reveal conspicuous sequence homogeneity, environmental isolates of B. cereus and B. thuringiensis exhibit extensive genetic diversity. Here we report the sequencing and comparative analysis of the genomes of two members of the B. cereus group, B. thuringiensis 97-27 subsp. konkukian serotype H34, isolated from a necrotic human wound, and B. cereus E33L, which was isolated from a swab of a zebra carcass in Namibia. These two strains, when analyzed by amplified fragment length polymorphism within a collection of over 300 of B. cereus, B. thuringiensis, and B. anthracis isolates, appear closely related to B. anthracis. The B. cereus E33L isolate appears to be the nearest relative to B. anthracis identified thus far. Whole-genome sequencing of B. thuringiensis 97-27and B. cereus E33L was undertaken to identify shared and unique genes among these isolates in comparison to the genomes of pathogenic strains B. anthracis Ames and B. cereus G9241 and nonpathogenic strains B. cereus ATCC 10987 and B. cereus ATCC 14579. Comparison of these genomes revealed differences in terms of virulence, metabolic competence, structural components, and regulatory mechanisms.

Use of a real-time PCR TaqMan assay for rapid identification and differentiation of Burkholderia pseudomallei and Burkholderia mallei.

A TaqMan allelic-discrimination assay designed around a synonymous single-nucleotide polymorphism was used to genotype Burkholderia pseudomallei and Burkholderia mallei isolates. The assay rapidly identifies and discriminates between these two highly pathogenic bacteria and does not cross-react with genetic near neighbors, such as Burkholderia thailandensis and Burkholderia cepacia.